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I.  rNTUOiniCTlON 


Missile  warhoads  arc  ly  (luite  different  in  dcsi)>n  and 

construction  than  warlioads  used  in  other  military  applications,  such  as 
artillery  shells.  This  difference  mainly  is  a consequence  of  the  low 
launch  loads  oxiicricnced  liy  missile  warhoads.  Warhead  design  In 
apjU  i cat  ions  where  high  launch  loads  are  experienced  (artillery  shells) 
arc  usually  of  the  natural  fragmentation  or  controlled  fragmentation 
type.  By  way  of  background,  natural  fragmentation  warheads  are  continuous 
metal  casing  warheads  with  no  preselected  control  sites  of  failure. 
Controlled  fragmentation  warheads,  on  the  other  hand,  are  continuous 
warheads  with  preselected  control  sites  or  lines  of  failure.  Those 
failure  sites  are  usually  arrived  at  mechanically,  mctallurgically  or 
explosively.  However,  missile  warheads  are  generally  discrete-fragment 
warlieads  and  not  continuous  metal-casing  systems.  Iliscrete-fragment 
warhoads  arc  preformed -fragment  warheads  in  which  the  fragments  of  the 
desired  shape  and  mass  are  placed  on  a liner  (metal  or  plastic)  which 
encases  the  high  explosive.  The  warhead  is  then  placed  Inside  a missile 
casing  (skin).  Those  proformed-fragments  are  attached  to  the  liner  and 
neighboring  fragments  by  an  epoxy  filler.  Due  to  the  very  low  launch 
load  environment  of  a missile,  this  evidently  is  sufficient  for  structural 
integrity. 

In  the  past  the  main  analytical  method  for  predicting  fragmentation 
warhead  performance*,  i.e.,  fragment  speed  and  direction,  has  been  the 
1 2 

Gurney  and  Taylor  formulae.  These  formulae  are  one-dimensional  by 
assumption  and  the  Gurney  formula  requires  an  empirically  determined 
constant.  This  constant,  called  the  Gurney  constant,  is  obtained 
experimental ly  via  warhead  tests  for  each  different  situation.  Due  to 
the  one-dimensional  nature  of  these  formulae,  their  utility  for  predicting 
warhead  performance  near  the  ends  is  questionable.  It  has  been 

3 

shown*  that  for  cases  when  the  flow  is  essentially  one-dimensional,  such 


*War>head  performance  as  used  in  this  report  will  mean  fragment  velocity 
distributio-n,  which,  it  can  be  shown,  yields  fragment  maos  distribution. 

Gurney,  /f.  W,,  "The  Initial  Velocities  of  Fragments  from  Bombs,  Shells 
and  Grenades",  BEL  Report  40b,  September  J04S,  AD^SOSlB, 

^‘Birkhoff,  G.,  MaaDougalt,  P.,  Pugh,  E,  M.  and  Taylor,  G.  I.,  Sir., 
"F,.vplo8ive8  with  Lined  Cavities",  <J.  Appl.  Phys.,  Vol.  19,  p.  663, 

June  1948. 

^Karpp,  H.  H.  and  Predebon,  W.  W.,  "Calculations  of  Fragment  Velocities 
from  Naturally  Fragmenting  Munitiona",  BEL  Memorandum  Report  No.  2609, 
July  1976.  (AD  #B007.377L) 
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as  lonn  iirtillory  projectiles  (L/D  > 2),  these  formulae  are  adequate 
for  predict ina  fraaraont  velocity.  However,  for  cases  where  the  flow  is 
two-dimensional,  such  as  near  the  warhead  ends  for  any  size  projectile, 
and  for  axi symmetric  warheads  with  length  to  diameter  ratios  less  than 
or  equal  to  two,  those  formulae  are  inadequate^.  Lastly,  these  formulae 
were  initially  derived  for  continuous  warheads  (i.o.,  natural  fragmenta- 
tion warhoadsj  and  for  one  specific  initiation  posture.  Consequently, 
for  discretc-frugment  warheads,  these  formulae  are  not  applicable  with- 
out significant  modification  and  experimental  recalibration. 

In  an  effort  to  find  a method  to  predict  warhead  performance  not 
only  over  its  midsoction  but  over  its  entire  length  including  the  ends, 
a time-dependent,  two-dimensional,  Lagrangian  computer  code  has  been 
applied  to  naturally  fragmenting  warheads.  In  Reference  3,  finite- 
difference  calculations  of  warheads  with  varying  explosive  fills  and 
casing  materials  wore  compared  with  experimental  results;  good  agreement 
was  demonstrated.  However,  the  computer  code  used  to  obtain  the  cal- 
culations in  Reference  3 assumes  a continuum.  Warheads  eventually 
fragment,  and  after  fragmentation  the  continuum  assumption  needs  to  be 
modified.  A method  to  model  the  effects  of  fragmentation  and  subsequent 
explosive  gas  leakage  was  presented  in  Reference  3 as  an  integral  part 
of  the  computational  scheme. 

In  this  report  the  earlier  work  by  Karpp  and  Predebon^  on  natural 
fragmenting  warhead  is  extended  to  missile  warheads,  i.e.,  discrete- 
fragment  warheads.  The  objective  of  the  work  reported  here  was  to 
develop  u warhead  modeling  capability  to  predict  the  performance  of 
missile  warheads  in  general,  as  well  as  warhead  designs  pertinent  to  the 
"Pocused-Blast-Fragment"  advanced  development  program,*  Computational 
modeling  was  emphasized  with  limited  experimental  verification  as  needed. 

In  tho  next  section  the  experimental  program  is  discussed  in  detail. 
In  Section  III  the  computational  modeling  is  explained.  Comparisons 
between  the  computational  and  experimental  results  are  discussed  in 
Section  IV.  The  conclusions  and  recommendations  follow  in  Section  V. 


*Thia  work  wae  aupported  hy  the  BRLj  Piaatinny  Areenal  and  AMC,  HQ. 
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II.  i:xpi:iumi;nta!.  pkocram 


I'lgurc  1 i.s  a schematic  drawint;  of  part  of  a tyi)ical  test  set-up 
used  in  the  warhead  test  firings  for  this  program.  .Shown  is  the  warhead 
on  an  angied*  wood  stand  which  rests  on  metal  plates  on  the  t loor  of  ttio 
seeond  level  of  tlie  test  facility.  The  two  metal  plates  form  a preset 
aperture  through  wliich  passes  the  column  of  fragments  to  bo  radiographed. 
The  lower  level  of  the  test  facility  contains  two  ortliogonal  banks  (each 
hank  has  four  flash  x-ruy  tubes)  which  provide  timc-scquentiul , ortho- 
gonal radiography. 

Uuo  to  the  funding  and  time  constraint.s  on  part  of  thi.s  program,  it 
becanio  obvious  at  tlie  very  beginning  that  cubical  fragments**  could  not 
be  placed  completely  around  each  tost  warhead  and  still  complete  the 
minimum  number  of  warhead  tests  planned,  Also,  from  a long  term  point 
of  view,  it  is  desirable  to  have  a more  economical  method  for  testing 
jn’oformed- fragment  warheads.  The  method  used  determined  the  minimum 
acimuthal-soctor  angle  containing  preformed- fragments  which  adequately 
represented  the  situation  where  the  entire  cylindrical  warhead  contained 
preformed- fragments.  All  the  warheads  tested  had  an  overall  cylindrical 
shape.  Figure  2 Illustrates  a tost  warhead  with  an  azimuthal  frogment- 
hcaving  sector  angle  of  22.5*  for  a typical  length  to  diameter  (1</D) 
ratio  of  2.0.  In  all  tests  whore  an  azimuthal  sector  of  the  warhead 
contained  the  preformod-fragmonts,  the  remainder  of  the  warhead  con- 
tained strips  (rods)  of  the  same  material.  The  width  of  the  strip.s  wa.s 
equal  to  the  width  of  the  fragments  (for  proper  radial  explosive  gas 
leakage),  and  the  strip  lengths  wore  equal  to  the  explosive  charge 
length. 

Table  I contains  all  the  characteristics  of  the  test  warheads  by 
their  respective  round  number.  As  indicated  in  Table  I,  the  warhead 
length  to  diameter  ratio  was  varied  from  1.0  to  2.0.  The  explosive 
charge  in  all  cases  was  Octol  (25 » TNT,  75‘i  IIMX).  The  fragments  wore 
steel  (ATSI  1018  or  1020)  cubes  and,  in  all  but  one  test,  wore  placed 
directly  on  the  bare  exj^losivc  charge.  This  typo  of  configuration  was 
adopted  in  order  to  determine  with  greater  accuracy  the  effect  which 
explosive  gas  leakage  between  fragments  has  on  the  fragment  velocity 
distribution,  particularly  for  the  computational  modeling  portion  of  the 
]irogram. 


Figures  3 and  4 illustrate  tyi^ical  radiographic  results  from  l./l)  - 
1.0  and  h/Il  = 2.0  warheads,  respect ivoly.  Both  tost  warheads  wore  initi- 
ated on  tlie  axis  at  the  loft  end.  Figures  3 and  4 show  the  longitudinal 
and  orthogonal  views  of  the  fragment  distributions  in  their  actual  spaclal 
positions  at  two  different  times.  In  each  of  the  tost  cases  the  fragments 
in  the  radiographed  column  were  premarked  for  identification  and  wore 
recovered  in  Colotex.  The  spacial  location  of  the  recovered  premtirked 

^The  angled  wood  stand  ie  need  to  rotate  the  main  fragment  distribution 
into  the  available  radiographia  view. 

**In  all  the  warhead  testa  the  fragment  shape  was  oubieal. 
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Figure  1.  Typical  Test  Configuration 
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TABLE  I-  aiARACTERTST ICS  AND  TT-ST  HATA 


i’rat;monts,  In  conjunction  with  the  t i me- sequent  ial , orthogonal  radiography 
ol‘  the  fragments,  permitted  dovclopnient  of  the  results  shown  in  f igures  5 
and  C for  the  two  tests  illustrated  in  figures  S and  4.  In  figures  5 and 
o,  each  of  the  fragments  in  flight  arc  labeled  according  to  their  original 
position  on  the  warhead  {shown  at  the  top  of  the  figures).  Two  different 
times  arc  shown  for  both  longitudinal  and  orthogonal  views. 

To  illustrate  rad  iogruiihical  ly  the  fragment  focusing  that  can  l)c 
achieved,  figure  7 shows  (for  two  different  times)  the  longitudinal  and 
orthogonal  views  of  an  L/D  = 2.0  warhead  which  was  Initiated  simultan- 
eously at  both  ends  on  the  axis  of  symmetry.  The  fragments  of  figure  7 
uro  labeled  in  Figure  8 according  to  their  original  position  cm  the 
warhead.  Again  two  different  times  arc  shown  for  the  longitudinol  and 
orthogonal  views. 

To  determine  the  accuracy  of  the  experimental  procedure,  a simple 
error  analysis  of  the  measurements  was  conducted  and  the  round  to  round 
variation  was  quantified.  Using  a conservative  approach  to  error 
analysis  wherein  only  maximum  difforonccs  are  considered,  for  two  time- 
sequential  observations  of  fragments  the  measurement  uncertainty  in  the 
fragment  speed  is  less  than  la  and  is  loss  than  * 0.004  for  the  tangent 
of  the  fragment  projection  angle.  The  round  to  round  variation  in  the 
warhead  tests  is  shown  in  Figure  9.  Plotted  in  Figure  9 is  the  fragment 
speed  and  projection  angle  distributions  vor.su.s  fragment  number  or 
relative  initial  axial  position  from  two  tests  of  nominally  identical 
(See  Table  I)  L/Il  = 2.0  cylindrical  warheads.  The  round  to  round 
variation  in  fragment  speed  Is  ± 3*  for  two  time-sequential  observations 
of  the  fragments  and  ± 13-o  for  single  time  observations  of  fragments 
(which  occurs  for  a few  fragments  at  both  ends).  The  overall  round  to 
round  variation  in  the  tangent  of  the  fragment  projection  angle  is 
•1  0.02.  Thoreforo,  it  appears  that  the  variations  from  one  round  to 
another  (nominally  the  sumo)  are  larger  than  the  measurement  uncertainty. 
Consequently,  it  was  inadvisable  to  attempt  to  improve  the  accuracy  of 
the  measurement  procedure. 

As  previously  mentioned,  the  cubical  fragments  were  placed  on  the 
warhead  in  varying  azimuthal  sector  angles,  with  strips  in  the  remaining 
azimuthal  angle,  to  determine  the  minimum  azimuthal  fragment -bearing 
sector  angle  required  to  approximate  the  situation  in  which  fragments 
covered  the  entire  round.  Plotted  in  Figures  10  and  11  are  the  fragment 
speed  and  projection  angle  distributions  for  varying  azimuthal  fragment- 
bearing  sector  angles  a for  L/0  = 1.0  and  2.0  warheads,  respectively. 

Since  the  percent  difference  between  the  a = 45®  and  a = 90®  cases  in 
figure  10  is  within  round  to  round  variation  shown  in  figure  9,  it  is 
felt  that  a 45“  azimuthal  fragment -bearing  sector  angle  is  adequate  for 
tho  L/D  = 1.0  case,  for  the  L/D  =2.0  case,  figure  11  indicates  a large 
difference  between  oi  = 22.5“  and  a = 45®,  90“  and  135®.  However  the 
difforoncos  are  much  smaller  between  the  a = 45®,  90®  and  135"  cases, 
although  the  fragment  speeds  for  a = 45°  tend  to  be  greater  than  or 
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figure  8.  Identification  of  Fragments  in  Figure  7 
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Ficuro  11.  Variation  of  "Warhead  Performance"  with  Aiimuthul 
Fragment -Bearing  Section  Anglo,  Wl)  *2.0 
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equal  to  u = 90"  or  a = 135°,  espcciully  In  areas  remote  from  otul  rure- 
faetlons.  Mindl'ul  of  the  latter  point  and  the  round  to  round  variation 
shown  in  Figure  9,  the  variation  in  fragment  speed  between  the  a - 90° 
and  135°  is  within  the  round  to  round  variation;  whereas  the  same  varia- 
tion lietween  the  a = 45°  and  90°  is  generally  not.  In  view  of  the  above 
it  is  felt  that  u = 90°  is  desirable  for  tlie  I./H  = 2.0  case. 

Returning  to  the  a = 22. E°  case  in  Figure  11,  due  to  the  unexpected 
largo  differences  between  it  and  a “ 45°,  90®  and  135°,  particularly 
near  the  ends  where  one  would  expect  the  fragment  speeds  to  converge  to 
a common  value  for  most  practical  it's,  it  is  questionable  whether  these 
large  diffoionces  are  duo  to  incroosed  confinement  for  a •=  22.5°  or 
perhaps  some  uncovered  error  in  the  data.  An  additional  teat  ut  a = 

22.5°  will  be  required  to  clarify  those  questions;  however,  it  wa.s 
judged  that  the  results  of  such  a test  would  not  significantly  change 
the  main  thrust  and  purpose  of  this  study  and  therefore  was  not  com- 
pleted ut  this  time. 

Tho  experimental  results  of  cases  in  which  computational  results 
wore  obtained  are  discussed  in  Section  IV.  There  arc  a few  test  cases 
for  which  tho  computations  were  not  performed  during  this  program  and 
those  experimental  results  will  be  presented  at  this  time. 

Tho  modeling  of  multi-fragment  layering  is  immensely  complicated 
duo  to  the  many  combinations  of  layer  patterns  that  can  bo  placed  on  the 
warhead.  In  order  to  bracket  the  effects  of  tho  simplest  multi- fragment 
layering  case,  l.e,,  two  layers  of  fragments,  two  double-layered  discrete- 
fragment  warheads  were  fired.  Tho  fragment  layers  were  aligned  with 
respect  to  ouch  other  for  tho  first  warhead;  the  other  warhead  hud  the 
fragments  bricked  (that  is,  the  fragments  on  the  top  layer  were  di.splaced 
by  a distance  equal  to  half  of  a cube  side  relative  to  tho  fragments  on 
the  bottom  layer).  In  Figure  12  are  plotted  the  results  of  tho  aligned 
and  bricked  double-layered  warheads.  In  order  to  compare  these  results 
to  tho  single  layered  test,  the  charge  mass  to  metal  muss  (C/M),  the 
typo  of  explosive,  explosive  initiation  and  explosive  dimensions  were 
kept  tho  same  us  tho  single  layered  1,/1)  2,0  case  (shown  in  Figure  16). 

F'urther,  tho  fragment  shape  and  fragment  material  wore  kept  tho  same. 

Those  constraints  thus  imposed  that  the  double- layered  fragment  musses 
be  1/8  the  single  layered  masses.  Due  to  the  resulting  lower  cube 
masses  (0.486  gm) , and  therefore  much  smaller  dimensions,  tho  previously 
u.sed  cube  numbering  system  was  not  effective  in  identifying  the  cubes 
with  respect  to  their  original  position  on  tho  warhead.  When  the  cubes 
were  recovered  they  had  generally  collided  with  ouch  other  during  flight 
and  tho  punched  numbers  on  each  fragment  were  unidontlfiablo.  Consequently 
the  fragment  speed  is  plotted  with  respect  to  its  projection  angle  in 
Figure  12  rather  than  its  initial  position  on  the  warhead. 
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Figure  12.  Double-Layered  Discrote-Frugmont  Warheads: 

Fragment  Speed  Versus  Projection  Angle 


In  I'liiurc  12  it  is  nuite  evident  that  there  exists  a significant 
variation  of  the  magnitude  of  tlic  velocity  at  any  projection  angle;  this 
oliscrvation  applies  both  to  the  aligned  and  bricked  warheads.  By  compari- 
son, if  the  data  from  the  single  fragment  layer  case  of  figure  16  were 
plotted  versus  projection  angle,  as  in  figure  12,  only  one  fragment  would 
appear  at  any  particular  projection  angle.  It  is  important  to  note,  for 
reason.s  previously  discussed,  it  was  very  difficult  to  identify  the  same 
fragment  in  two  time-sequential,  radiographic  views  for  many  of  the  frag- 
ments in  the  double-layered  cases.  Consequently,  the  only  fragments  that 
arc  plotted  are  those  in  which  this  identification  was  confidently  made. 
Therefore,  although  the  single  fragment  layer  case  cannot  be  compared 
with  the  double- layered  case  over  the  entire  projection  angle  distribution, 
meaningful  comparisons  can  be  made  over  the  projection  angle  range  for 
which  data  wore  obtained.  In  the  two  examples  depicted  in  Figure  12,  the 
spread  in  the  speeds  ranges  from  0.15  to  0.40  mm/psec  over  the  projection 
angle  distribution.  The  mo.st  probable  contributing  cause  of  this  spread 
in  speed  at  n given  projection  angle  is  slippage  between  the  layer  frag- 
ments. Thus,  instead  of  fragments  leaving  the  warhead  at  a given  angle 
having  one  particular  velocity,  collective  effects  result  in  a number  of 
fragments  with  varying  velocities  at  the  specified  angle. 

It  is  lastly  noted  that  the  plotted  speeds  in  Figure  12  have  not 
been  corrected  for  drug,  which  from  the  data  was  estimated  to  be  causing 
a 4'd-7"o  decrease  in  fragment  speeds.  This  decrease  wn.s  obtained  by 
comparing  single  and  double  flash  calculations  of  Identical  cubes.  This 
calculation  was  also  completed  in  the  single  layered  cases;  however,  in 
these  cases  the  drag  was  always  negligible.  Because  of  the  lower  cube 
masses  (0.486  gmj  in  the  double-layered  warheads,  drag  is  no  longer 
negligible. 

Sho\m  in  Figure  13  are  the  experimental  results  for  single-layered 
fragments  on  an  aluminum  liner  with  aluminum  endplates,  L/D  = 2.0  cylin- 
drical warhead.  As  shown,  initiation  was  on  the  axis  on  the  left  side  of 
the  warhead.  All  the  necessary  data  are  given  in  Table  I.  This  is  the 
only  test  case  reported  with  a liner  and  endplates.  These  results  will 
be  discussed  in  some  detail  in  Section  IV  along  w'-h  a similar  test  case 
without  a liner  and  endplates. 


III.  COMPUTATIONAL  MODELING 

The  computer  code  that  is  being  used  for  the  computation  modeling  is 
a modified  version  of  the  HliMP'*  computer  code.  The  UF.MP  code,  which  was 
originally  developed  at  the  Lawrence  Livermore  Laboratories,  California, 
is  a time -dependent , two-dimensional,  Lagrangian  finite-difference  code. 

''WilkinOj  M.  L.j  "Calaulationa  of  Etaatia-Plaetia  Flow^"  Methods  of 
Computational  Phyoiae,  Vot.  3,  edited  by  Aider,  B. , Fernbaak,  S.  and 
Rotenburg,  M.,  Aaademia  Preao,  NY,  1964, 
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If  the  stundard  form  of  the  IIHMP  code  is  used  to  predict  the  performance 
of  discrotc-friigment  warheads  without  including  the  effects  of  fragment 
separation  and  subsequent  explosive  gas  leakage  between  fragments,  the 
dashed  line  results  shown  in  Figure  14  are  obtained.  Shown  in  Figure  14 
is  a comparison  between  these  HEMP  calculations  and  experimental  data  of 
fragment  speed  and  projection  angle  distributions  for  an  L/D  =1.0  cylin- 
drical warhead.  In  the  HEMP  calculations  the  fragments  are  modeled  as  a 
fluid,  i.o.(  the  fragments  cannot  accommodate  a tensile  force  but  they 
can  be  recompressod.  Ihis  would  be  a logical  approach  with  a continuum 
code  since  the  fragments  are  not  a continuum  but  discrete  particles  which 
will  separate  with  the  application  of  a small  tensile  force.  It  is 
obvious  that  this  modeling  Is  superior  in  predictive  ability  than  the 
prediction  of  velocity  and  angle  by  Gurney  and  Taylor,  also  shown  in 
Figure  14.  However,  the  agreement  between  the  HEMP  computations  and 
experiment  is  not  as  good  as  desired.  Because  of  this  and  other  objec- 
tions®, modeling  of  the  effects  of  fragment  separation  and  explosive  gas 
leakage  between  fragments  was  initiated  to  improve  the  overall  predictive 
capability  of  the  computations. 

The  discrete-fragment  warhead  (missile  warhead)  model  presented 
Involves  the  following  three  essential  features:  A)  modeling  of  the 
discrete-fragments,  B)  modeling  of  fragment  separation  and  subsequent 
explosive  gas  leakage  and,  C)  criteria  for  fragment  separation  and 
method  to  determine  this  criteria. 

A.  Modeling  of  Discrete-Fragments 

When  discrete-fragments  are  placed  on  a cylindrical  body,  whether  it 
bo  the  bare  explosive  or  a liner  containing  the  explosive,  (assuming 
conventional  fragment  shapes)  air  gaps  will  occur  between  the  fragments 
in  the  circumferential  direction.  Consequently,  the  discrete-fragment 
casing  forms  a composite  of  metal  and  air.  In  reality,  this  is  a three- 
dimensional  system.  However,  the  approach  that  was  adopted  with  the  two- 
dimensional  modeling  is  similar  to  the  methods  used  in  the  modeling  of 
composite  materials.  The  discrete  metal-air  composite  casing  (metal 
preformed-fragments)  is  modeled  as  a homogeneous  continuous  material  of 
the  same  metal  material  as  the  experimental  fragments,  but  with  the 
requirement  that  the  total  mass  of  the  new  continuous  motal  casing  equal 
the  total  miss  of  all  the  fragments  on  the  experimental  warhead.  For 
example,  in  the  case  of  cubical  fragments,  the  resulting  continuous  metal 
casing  wall  thickness  in  the  model  will  be  less  than  one  of  the  sides  of 


^Karppj  E.  E.  and  PredeboKj,  W.  W.^  "Caloulationa  of  Fragment  VeloaitieB 
from  Fragmentation  Munitionej  " First  International  Sympoaium  on  Bailie- 
tioB,  Orlando,  Florida,  IB-IS  November  1974,  Section  IV,  pp.  146-176, 
Prooeedinga  Ftibliehed  by  American  Defense  Preparedneee  Aeaoaiation,  Union 
Trust  Building,  Vaahington,  DC. 
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an  actual  cubical  fiaHmcat . I'or  a cylindricul  warhead,  the  modeled 
casiiiji  outside  diameter  Tor  a fjlven  explosive  diameter  is  obtained  from 

M,|,  = 1'^.  1.^  “ ‘’lll/^  * * 

1 

where  M.j,  is  the  measured  total  fragment  mass  (plus  liner  mass  if  the 

warhead  has  a linerl,  iiij  are  the  mousured  individual  fragment  masses,  p^, 

is  the  fragment  density,  1.^^  is  the  total  length  of  the  fragments  on  the 

warhead,  llm,  is  the  measured  diameter  of  the  explosive  and  llj^  Is  the 

outside  diameter  of  tlie  continuous  mctul  casing  to  bo  computed.  If  one 
ignores  the  ulr  gap  between  the  fragments  and  uses  the  fragment  cube 
dimension  us  the  continuous  casing  wall  thickness,  thou  for  the  cases 
reported  in  Sections  11  and  IV,  the  inodolod  total  mass  will  bo  us  much 
a.s  greater  than  the  actual  total  muss. 

B.  i'ragment  Separation  and  lixploslvc  Gas  Leakage 

The  general  Ideas  of  the  model  for  fragment  separation  and  explosive 
gas  leakage  will  he  presented  hero;  however,  the  specific  details  and 
oiiuatlons  have  been  presented  in  an  earlier  report®. 

The  llliMP  code  models  continuum  behavior  In  general.  However,  in 
reality,  when  fragment  separation  occurs  a continuous  casing  circumfer- 
ence cca.scs  to  increase.  But  the  standard  form  of  the  code  causes  the 
casing  to  remain  continuous.  This  treatment  tends  to  yield  a higher 
radial  acceleration  and  thus  a higher  final  fragment  velocity.  This  cun 
bo  corrected  by  computing  the  true  force  acting  on  an  element  of  the 
casing.  The  true  force  is  computed  In  the  code  by  reducing  the  pressure 
acting  on  the  interface  of  the  casing  by  the  factor  where  is 

the  radius  at  fragment  separation  and  R is  tho  current  radius.  At 
present,  must  be  determined  by  experiment.  After  fragment  separation, 

i.e.,  U > a void  appears  between  fragments  and  gas  leakage  occurs. 

i'bis  effect  is  modeled  In  the  code  by  calculating  tho  rate  of  efflux 
from  an  ideal  nozzle  passing  its  maximum  rate  of  flow.  The  rate  of 
efflux  is  computed  in  the  code  for  every  cycle,  and  the  mass  leaked  for 
a single  cycle  is  the  product  of  the  muss  rate  and  the  cycle  time 
Increment.  In  this  manner,  the  model  accounts  for  gas  leakage  between 
fragments  after  casing  breakup.  Hus  leakage  reduces  tho  pressure  accord- 
ingly. Lastly,  to  simulate  a loss  in  circumferential  stress  upon  casing 
breakup,  the  yield  strength  of  the  fragment  material  is  set  equal  to 
zero  vjhen  R > 

C . Criterion  for  I'ragment  Separation 

Currently,  the  criterion  for  fragment  separation,  R^^,  is  an  empiri- 
cally determined  quantity.  As  discussed  in  Section  II  the  fragments 

SI 


arc  all  prcmnrkL'd  and  recovered  in  a soft  recovery  medium  (Celotcx) 
after  firing.  The  recovered  fragments  arc  weighed  and  measured  for 
their  dimensions.  Obtained  from  these  moasuremonts  is  the  fragment 
plastic  deformation  due  to  the  explosive  loading  (or  final  dimensions^. 
Since  its  original  dimensions  are  known  a straight  forward  calculation 
determines  the  corresponding  radius  of  the  warhead  at  which  the  fragments 
separated.  (The  elastic  displacements  are  neglected  in  this  analysis), 
for  the  case  of  cubical  fragments  on  bare  explosive  for  L/D  * 1.0  and 
2.0  cases,  the  expansion  ratio  at  fragment  separation,  i.e.,  where 

is  the  inside  radius  of  the  casing,  was  found  to  bo  1.2.  The  recovered 

fragment  measurements  and  expansion  ratios  are  tabulated  In  Table  1. 

This  criterion  for  fragment  separation  is  empirical;  the  criterion  will 
not  be  applicable  in  general.  Thus,  the  next  refinement  in  this  analysis 
is  to  relate  this  criterion  to  a fundamental  material  property. 


IV.  COMPARISON  BETWIiKN  COMPUTATIONAL  AND  EXPERIMENTAL  RESULTS 

In  Figures  15-17  are  shown  comparisons  between  the  calculations 
with  the  discrete-fragment  model  (labeled  Elastic-Plastic  with  Gas 
Leakage  Model)  and  experimental  ro.sult3.  In  Figure  IS  and  16  the  calcu- 
lated fragment  speed  and  projection  angle  distributions  are  compared 
with  the  experimental  data  for  L/D  - 1.0  and  2,0  cylindrical  warheads, 
respectively.  Data  from  two  test  firings  are  plotted  in  each  figure. 

In  both  cases  the  agreement  between  computation  and  experiment  is  extremely 
good.  Also  shown  in  Figures  IS  and  16  are  the  calculations  without  the 
discrete -fragment  model,  namely,  a fluid  model,  and  the  Gurney  velocity 
and  Taylor  angle  predictions.  Shown  in  Figure  17  is  a comparison  between 
the  computations  and  experimental  data  for  the  same  case  as  shown  in 
Figure  16  except  with  dual-end  axial  initiation.  Again  the  agreement 
between  computation  and  exporimont  is  very  good.  Note  that  the  computa- 
tions model  the  region  where  the  two  detonation  waves  collide  (higher 
pressure  region)  at  0.5  relative  initial  axial  position. 

Returning  to  the  tost  results  shown  in  Figure  13,  note  that  the 
warhead  configuration  in  Figure  13  is  similar  to  the  configuration  in 
Figure  16.  Both  are  L/D  = 2.0  cylindrical  warheads  with  the  same  explo- 
sive material,  fragment  material  and  size,  and  with  essentially  the  same 
explosive  diameter.  However,  the  test  warhead  in  Figure  13  has  a liner 
between  the  explosive  and  fragments,  and  ondplatos.  The  reasons  given 
for  placing  the  liner  between  the  fragments  and  the  explosive  are  often 
varied.  It  was  the  intent  of  this  initial  experiment,  in  conjunction 
with  computational  modeling  (yet  to  be  performed),  to  help  e.\plain  the 
effect  of  a liner  on  warhead  performance.  Duo  to  the  additional  width 
of  the  A1  liner,  the  inside  fragment  radius  increased,  which  resulted  in 
more  fragments  on  the  warhead,  and  thus  a larger  total  metal  weight,  and 
lower  C/M.  These  data  are  given  in  Table  1 by  Rd.  No.  9457.  Consequently, 
the  results  of  Figures  13  and  16  cannot  bo  directly  compared.  However, 
it  is  the  authors'  judgment  at  this  time,  based  on  these  data  and  other 
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l'i>!ure  15.  Comparison  of  Fluid  Model  Computations  and  Dis- 

croto-Fragmont  Model  Computations  with  Experimental 
Data,  1./D  “1.0  Cylindrical  Preformed -Fragment 
Warhead 
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Figure  16.  Comparison  of  Fluid  Model  Computations  and  Di-s- 

creto- Fragment  Model  Computations  with  Experimental 
Data,  L/D  ■2.0  Cylindrical  Preformed- Fragment 
Warhead 
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calculut ions,  that  the  effect  of  u liner  Is  to  contain  the  explosive 
Hiisos  for  a lonficr  period  of  time.  This  in  turn  causes  leakage  of  explo- 
sive gases  between  fragments  to  occur  at  a larger  expansion  ratio,  result 
Ing  in  u higher  final  fragment  speed,  as  compared  to  the  case  without  a 
liner  (assuming  all  other  quantities  being  equal).  This  can  be  illus- 
trated from  a plot^  of  calculated  fragment  speed  versus  time  (from  dotona 
tion)  for  those  cases.  1-rom  a plot  of  this  typo  one  observes  that  frag- 
ment separation  and  leakage  occurs  during  the  initial,  approximately 
exponential,  acceleration  rise.  By  placing  a liner  between  the  fragments 
and  explosive,  one  can  prolong  fragment  separation  and  subsequent  leakage 
until  the  fragment  speed  versus  time  profile  flattens,  which  then  results 
in  a higher  final  fragment  velocity,  from  those  and  other  computations, 
it  is  estimated  that  up  to  approximately  a 10%  increase  in  final  fragment 
velocity  can  bo  obtained  with  a liner  between  the  fragments  and  explo- 
sive versus  the  case  without  a liner. 

Having  demonstrated  extremely  good  agreement  between  the  present 
calculations  and  experimental  data  for  both  L/D  “ 1.0  and  2.0  cylindrical 
warheads  with  single-end  and  dual-end  axial  initiation,  the  computations 
wore  then  confidently  utilized,  without  experimental  verification,  to 
study  the  effect  of  initiation  posture  on  warhead  performance.  Shown 
in  Figure  18  are  the  nine  initiation  postures  that  were  computed,  The 
cylindrical  warheads  that  are  drawn  in  Figure  18  are  all  with  a length 
(L)  to  diameter  (0)  ratio  of  two  and  initiated  as  shown.  The  first  case 
of  dual  end  axial  booster  initiation,  l.e.,  initiation  with  two  19.05  mm 
diameter  explosive  boosters,  is  shown  in  Figure  17,  The  computational 
results  of  the  remaining  eight  cases  are  shown  respectively  In  Figures 
19-26.  Shown  in  Figure  17  and  Figures  19-26  is  the  effect  of  a large 
variety  of  initiation  postures  on  warhead  performance.  The  general  idea 
of  this  set  of  computations  is,  on  the  one  hand,  to  demonstrate  the 
versatility  and  generality  of  this  computational  tool,  and  on  the  other 
hand  to  aid  in  the  understanding  of  the  effect  of  initiation  posture  on 
warhead  performance.  Furthermore,  should  the  initiation  posture  of 
interest  lie  somewhere  between  the.se  cases,  a simple  interpolation  of 
the  already  computed  cases  is  all  that  is  required  for  a first  approxima- 
tion of  the  warhead  performance.  Similarly,  since  the  computational 
results  are  general,  i.e.,  scalable,  then  if  the  geometry  of  interest  is 
not  procl.sely  the  same  as  the  cases  calculated,  these  results  can  bo 
scaled  up  or  down  assuming  proper  scaling  procedures  are  followed. 

Lastly,  it  is  worth  while  to  note  that  axial  lino  initiation.  Figure  26, 
Is  the  one  case  of  the  initiation  variations  studied  in  which  the  frug- 
inont  projection  angle  and  speed  distribution  are  significantly  different 
from  the  dual  end  axial  initiation  case,  Figure  17. 

Shown  in  Figures  27  through  30  are  the  computational  results  for 
two  different  internal  cavity  sizes  and  materials  for  the  same  general 
L/D  2.0  cylindrical  geometry.  Shown  in  Figure  27  are  the  results  for 
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the  case  with  an  aluminum  cavity  liner  material  with  the  dimensions  as 
indicated.  Comparing  the  results  shown  in  Pigure  27  with  those  in 
Figure  17  it  can  bo  soon  that  although  the  projection  angle  distribution 
is  only  slightly  altered,  the  fragment  speed  is  significantly  reduced. 

The  computational  results  with  a tungsten  alloy  (p  = 16890  KG/M^) 
cavity  liner  material  for  the  same  cylindrical  geometry  and  cavity  shape 
given  in  Figure  27  are  shown  in  Figure  28.  As  expected,  due  to  the 
higher  density  cavity  liner  material  the  fragment  speed  distribution  is 
somewhat  higher  than  with  aluminum;  however,  the  fragment  projection 
angle  distribution  is  not  noticeably  altered.  It  is  Important  to  note 
that  these  results  are  for  one  cavity  shape  and  therefore  are  not  neces- 
sarily true  for  all  cavity  shapes  and  sizes,  In  Figures  29  and  30  are 
shown  the  computational  results  for  a slightly  different  cavity  size 
where  the  cavity  is  modeled  either  as  a vacuum  (Figure  29)  or  as  a 
cavity  whose  boundaries  are  rigid  (Figure  30)  for  all  time.  The  purpose 
of  these  two  computations  is  to  provide  a qualitative  lower  (Figure  29) 
and  upper  (Figure  30)  bound  on  fragment  speed  and  projection  angle 
distributions  that  can  be  expected  for  the  cylindrical  geometry  and 
cavity  shape  shown. 


V.  CONCLUSIONS  AND  RBCOMMENDATIONwS 

A di^cre^r-fragment  warhead  (missile  warhead)  computational  modeling 
capability  has  been  developed  which  is  capable  of  modeling  with  a satis- 
factory degree  of  confidence  the  performance  of  a large  variety  of  axi- 
symmetric  warhead  configurations,  The  computational  warhead  model 
Includes  the  ability  to  model  the  effects  of  varying  L/D,  fragment 
material  and  size,  explosive  material,  initiation  posture,  and  internol 
cavity  material  and  shape.  It  is  felt  that  the  modeling  capability  can 
be  extended  to  include  two  layers  of  fragments  on  the  warhead,  a liner 
between  the  fragments  and  the  explo.sive,  and  end  confinement.  The 
experiments  wore  completed  for  these  cases  and  are  discussed  in  Section 
11. 


Future  efforts  should  include  additional  experiments  and  computa- 
tional modeling  with  different  linor  materials,  shapes  and  sizes,  and 
end  confinement  to  understand  fully  the  Influence  of  these  variables  on 
warhead  performance.  Also,  future  work  should  include  generalization  of 
the  empirical  criterion  for  fragment  separation  to  a more  fundamental 
criterion.  For  example,  material  properties,  such  as  the  total  plastic 
strain  to  failure,  might  be  used  to  predict  fragment  separation. 

It  is  also  believed  that  the  modeling  capability  can  successfully 
be  applied  to  certain  aspect.s  of  eccentric  warheads  (a  three  dimensional 
problem).  However,  while  some  degree  of  success  may  be  had,  in  order  to 
model  eccentric  warheads  properly,  a throe  dimensional  computer  code  is 
needed. 


SO 
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L!ST  OF  SYMBOLS 


Nominal  charge  (explosive)  mass  to  metal  mass  (fragment  t 
liner)  ratio 

Explosive  diameter 

Calculated  outside  diameter  of  continuous  metal  casing  model, 
from  Mj,  ■ Pj  {ir/4  (Dj^j*  D^^g*)  } 

Side  of  a cubical  fragment 

Explosive  Length 

Length  of  fragments  on  warhead  (equal  to  L||g  In  all  cases) 

Mass  of  a cubical  fragment  (also  equivalent  to  m^) 

Mass  of  the  explosive  charge 

Mass  of  a rod 

Individual  fragment  mass 

Maas  of  a recovered  cubical  fragment 

Total  mass  of  fragments  and  liner  on  warhead 

Number  of  cubical  fragments  per  column  (along  warhead  length) 

Number  of  columns  of  cubical  fragments 

Number  of  rods  per  warhead 

Current  warhead  radius  in  computations 

Radius  of  warhead  at  fragment  separation  in  computations 

Expansion  ratio  at  fragment  separation  or  breakage  ratio; 
equivalent  to 

Initial  inside  radius  of  casing  (metal  fragments) 

Recovered  cube  dimension  in  the  circumferential  (hoop)  direction. 
Azimuthal  fragment-bearing  sector  angle 
Fragment  mass  density 


DISTRIBUTION  LIST 


No.  of-  No.  of 

Copica  Ortiuni  zut  ion  Copies  OrRoni  zatlon 


12  Commiindcr 

Defense  Documentation  Center 
ATTN:  DDC-TCA 
Cameron  Station 
Alexandria,  VA  22314 

1  Director 

Defense  Advanced  Research 
Projects  Agency 
ATTN:  Tech  Info 
1400  Wilson  Boulevard 
Arlington,  VA  22209 

1 Commander 

US  Army  Materiel  Development 
and  Readiness  Command 
ATTN:  DRCDMA-ST 
SOOl  Blsenhower  Avenue 
Alexandria,  VA  22333 

1 Commander 

US  Army  Materiel  Development 
and  Readiness  Conunand 
ATTN : DRCDL 
SOOl  Eisenhower  Avenue 
Alexandria,  VA  22333 

1 Commander 

US  Army  Aviation  Research 
and  Development  Command 
ATTN:  DRSAV-E 
12th  and  Spruce  Streets 
St,  Louis,  MO  63166 

1 Director 

US  Army  Air  Mobility  Research 
and  Development  Laboratory 
Ames  Research  Center 
Moffett  Field,  CA  94035 

2 Commander 

US  Army  Electronics  Command 
ATTN:  DRSEL-RD 

DRSEL-HL-CT,  S. Crossman 
Fort  Monmouth,  NJ  07703 


3  Commander 

US  Army  Missile  Research 
and  Development  Command 
ATTN:  DRDMI-R 
DRDMI-RPS 

DRDMl-RDP,  II. W.  Burnum 
Redstone  Arsenal , AL  35809 

1 Commander 

US  Army  Tank  Automotive 

Research  5 Development  Command 
ATTN:  DRDTA-RWL 
Warren,  MI  48090 

2 Commander 

US  Army  Mobility  Equipment 
Research  5 Development  Command 
ATTN:  Tech  Docu  Cen,  Bldg.  .315 
DRSME-RZT 

Fort  Boivoir,  VA  22060 
1 Commander 

US  Army  Armament  Materiel 
Readiness  Command 
Rock  Lsland,  IL  61202 

1 Commander 

US  Army  Armament  Research 
and  Development  Command 
ATTN:  DRSAR-CG 
Dover,  NJ  07801 

5 Commander 

US  Army  Armament  Research 
and  Development  Command 
AITN:  DRDAR-LC,  Dr.  J.  Frasier 
DRDAR-LCA, 

Dr.  E.  G.  Sharkoff 
Mr.  L.  Rosondorf 
Mr.  G.  Randers-Pehrson 
Mr.  J.  Pearson 
Dover,  NJ  07001 


53 


DiSTIUBUTION  LIST 


No . o t’ 

C.op  1 es  OrKitnizittion 


No,  of 

Copies  1 zution 


5 CommHiider 

US  A»’my  Armament  Itcsoarch 
and  Dovolopment  Command 
AITN:  DROAR-LCli,  Dr.  R. Walker 
DRDAR-TS-S 

DRDAR-LC,  Mr.  R.  Nixon 
Dover,  NJ  07801 

1 Commander 

US  Army  Hi>’ * v Diamond  Luba 
ATTN:  DRXDO-Tl 
2800  Powder  Mill  Road 
Adelphi,  MD  2078,' 

2 Commander 

US  Army  Materials  and 
Mechanica  Research  Center 
ATTN:  J.  Mescal 1 
D.  Roy lance 
Watertown,  MA  02172 

1 Commander 

US  Army  Natick  Research 
and  Development  Comm“nd 
ATTN:  DRXRH,  Dr.  D,  Sioling 
Natick,  MA  01762 

1 Director 

US  Army  TRADOC  Systems 
Analysis  Activity 
ATTN:  ATAA-SL  (Tech  Lib) 

White  Sands  Missile  Range 
NM  88002 

1 Deputy  Assistant  Secretary 
of  the  Army  (R8D) 

Department  of  the  Army 
Washington,  DC  20310 

1 HQDA  (DAMA-eSM) 

Washington,  DC  20310 

1 HQDA  (DAMA-ARP) 

Washington,  DC  20310 


1 HQDA  (DAMA-MS) 

Washington,  DC  20.310 

1 Commander 

US  Army  Research  Office 
ATl’N:  Dr.  li.  Saibol 
P.  0.  Box  12211 
Ro.search  Triangle  Park 
NC  27709 

1 Commander 

US  Naval  Air  Systems  Command 
ATPN:  AIR- 604 
Washington,  DC  20360 

1 Commander 

US  Naval  Ordnance  Systems  Cmd 
ATTN;  ORD-9132 
Washington,  DC  20360 

1 Chief  of  Naval  Research 

ATTN:  Code  ONR  439,  N.Perrono 
Washington,  DC  20360 

1 Commander 

US  Naval  Air  Development  Ctr 
Johnsvillo 

Warminster,  PA  18974 

1 Commander 

US  Naval  Missile  Center 
Point  Mugu,  CA  93041 

2 Commander 

David  W.  Taylor  Naval  Ship 
Research  8 Development  Center 
ATTN:  D.  R.  Garrison 
A.  Wilner 

Bothesda,  MD  20034 

3 Commander 

US  Naval  Surface  Weapons  Center 
ATTN:  Code  TBB,  D.W. Colbert son 
Mr , L , Hock 

Code  TX,  Dr,  W.G.  Sopor 
Dahlgron,  VA  22448 


54 


No.  of 
Con  i c .s 


DISTRIBUTION  LIST 


. 


1 1 

ii 


Oraan i zation 

2 Contmamler 

US  Naval  Surface  Weapons  Center 
ATTN:  WR-13 

Tech  Library 

Silver  Spring,  MD  20910 

6 Commander 

US  Naval  Weapons  Center 
ATIN;  J.  Pearson,  Code  383 

J.S.  Rinehart,  Code  383 
M.  Backman,  Code  3831 
C.F.  Austin,  Code  3833 
R.G.S.  Sewell,  Code  3835 
J.  Wooks,  Code  3266 
China  Lake,  CA  93555 

2 Commander 

US  Naval  Research  Laboratory 
ATTN ; Tech  Lib 
J . Baker 

Washington,  DC  20375 

1 Superintendent 

US  Naval  Postgraduate  School 
ATTN;  Dir  of  Lib 
Monterey,  CA  93940 

1 AFATL  [Mr,  Leonard  T.  Wilson) 
Eglin  AFB,  FL  32542 

3 AFATL  (DLRD,  K.  McArdle; 

DLRV,  Cl.  Crews;  DLYD) 

Eglin  AFB,  FL  32542 

1 ADTC/DLJW  (CPT  D.  Matuska) 
Eglin  AFB,  FL  32542 

] AFATL  (DLDL,  MAJ  J.E.  Morgan) 
Eglin  AFB,  FL  32542 


Organization 

2 AFML 

Wright-Patterson  AFB,  OH  454.33 
1 AFFDL  (FDT) 

Wright-Patterson  AFB,  OH  45433 

3 ASD  (YH/EX,  John  Riovley; 

XRHD,  Gerald  Bennett; 

BNYS,  Matt  Kolleck) 
Wright-Patterson  AFB,  OH  45433 

1 Headquarters 

National  Aeronautics  and 
Space  Administration 
Washington,  DC  20546 

2 Director 

National  Aeronautics  and 
Space  Administration 
John  F.  Kennedy  Space  Center 
ATTN:  KN-ES-3 
KN-ES-34 

Kennedy  Space  Center,  PL  32899 
2 Director 

National  Aeronautics  and 
Space  Administration 
Langley  Research  Center 
ATTN:  R.  J.  Hayduk 
Tech  Lib 

Hampton,  VA  23665 

1 US  Energy  Research  and 

Development  Administration 
Los  Alamos  Scientific  Lab 
ATTN:  Mr.  R.  Craig,  MS-960 
Dr.  L.  Mantel,  WX-2 
Dr.  Robert  Karpp 
Los  Alamos,  NM  87544 


No.  of 
Copies 


2 AFWL  (WLL) 

Kirtland  AFB,  NM  87117 

1 ASD  (ASBRS) 

Wright-Patterson  AFB,  OH  45435 

\ 

\ 

55 


\i 


